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SUMMARY 

A theoretical investigation of the velocity profiles for laminar 
mixing of a high-velocity stream with a region of fluid at rest 
has been made assuming that the Prandtl number is unity. 
A method which involves only quadratures is presented for 
calculating the velocity profile in the mixing layer for an arbitrary 
■value of the free-stream Mach number. 

Detailed velocity profiles have been calculated for free-stream 
Mach numbers of 0. 1, 2, 8, and 5. For each Mach number, 
velocity profiles are presented for both a linear and a 0.76- 
poiver variation of viscosity with absolute temperature. The 
calculations for a linear variation are much simpler than those 
for a 0.76-power variation. It is shown that by selecting the 
constant of proportionality in the linear approximation such 
that it gives the correct value for the viscosity in the high- 
temperature part of the mixing layer, the resulting velocity 
profiles are in excellent agreement with those calculated by a 
0.76-power variation. 
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INTRODUCTION , .* 

The relocity profile _for turbqlenj mixing at constant 
pressure of an incompressible' streaffi "with a dead-air region 
has been calculated by several investigators, principally 
Tollmien (reference 1). These calculations agree well with 
the available experimental data, although the conventional 
assumptions regarding the mixing length of a turbulent 
flow have since been shown by experiments to be incorrect 
(reference 2) . The many difficulties encountered in making 
precise turbulent- mixin g calculations are, of course, a 
consequence of the extremely complicated mechanism 
govering all turbulent flows. In contradistinction to the 
case of turbulent mixing, the mechanism involved in laminar 
mixing is relatively simple, and the mathematical relation 
between stresses and velocity gradients for laminar flow 
is well known. The velocity profiles for l amin ar mixing , 
however, apparently have not as yet been calculated even 
for the case of incompressible flow. It is the purpose of the 
present paper to calculate the velocity profiles for laminar 
mixing (starting with zero boundary-layer thickness) of 
an air stream of arbitrary temperature with a dead-air 
region also of arbitrary temperature. In cases where a 
laminar boundary layer of appreciable thickness exists at 
the point where mixing begins, the results given herein are 
not directly applicable in the initial part of the mixing 
region. For such cases, it is necessary to make some 
supplementary approximation in order to apply the results. 

Since’ the practical applications of laminar-mixing phe- 
nomenon 'usually involve the flow of a gas. the present 


analysis includes the effects of compressibility. Examples 
of typical flows wherein laminar mixing occurs can be found 
in the flow of small-scale jet pumps, in the flow behind the 
intersection of shock waves of unequal strength, and in 
the flow immediately behind the base of a body which has 
a laminar boundary layer. 1 

SYMBOLS AND NOTATIONS 

C constant of proportionality between viscosity anrl 
temperature 

c v specific heat at constant pressure 

1‘ coefficient of heat conduction 

L characteristic length 

M Mach number 

p static pressure 

Pr Prandtl number 

S Sutherlands constant, approximately 216° F for air 
T absolute temperature 

U<o free-stream velocity 

u, v velocity components in x, y directions, respectively 
x, y Cartesian coordinates 

5 thickness of mixing layer, taken between points where 
the velocity is 0.01 and 0.99 of the free-stream 
velocity 

co exponent of viscosity variation with temperature 
p mass density 

fi coefficient of viscosity 

v kinematic coefficient of viscosity 

f stream function 

£ dimensionless independent variable 

SUBSCRIPTS AND SUPERSCRIPTS 

* dimensionless variables as defined in equation (14) 

00 free-stream conditions 

0 stagnation conditions of the free stream 

d conditions in the dead-air region 

BASIC EQUATIONS AND ASSUMPTIONS 
BASIC EQUATION'S 

A schematic illustration of the flow under consideration is 
shown in figure 1. In order to make the l amin ar -mixing 
process amenable to calculation, the usual assumptions are 
made that the layer affected by viscosity is thin and has zero 

1 The present analysis was undertaken as part of an investigation of this latter problem, 
and originally appeared as Appendix B of a thesis “Base Pressure at Supersonic Velocities,” 
submitted to the California Institute of Technology, June 1948. The results of some supple- 
mentary computati* -ns not given in the thesis have been added for sake of completeness in 
the present report. 
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pressure gradient. Under these conditions the formal pro- 
cedure for estimating the order of magnitude of the various 
terms in the complete Navier-Stokes equations for viscous' 
compressible flow can be carried through in precisely the 
same manner as is done in the classical (Prandtl) treatment of 
laminar boundary-layer flows. In so doing, the dynamic 
equation for the ar direction reduces to the familiar boundary- 
layer momentum equation 

dtt , bu d / bu\ m 

pu bx +pv by~ by by) 

while the dynamic equation for the y direction reduces to 
zero on both sides. In passing from the Navier-Stokes equa- 
tion to (1) it is to be noted that the usual boundary-layer 
assumption 5/ce<< 1 is violated in the immediate vicinity 
of point 0 (fig. 1) just as in the case of boundary-layer flow 
over a plate. 

By employing the same considerations on order of magni- 
tudes as were used for the complete Navier-Stokes equations, 
the complete differential equation representing the balance 
of energy in viscous compressible flow reduces to 

b(e„T) , b(c P T) b bT\ , fbu\ ^ 

(2) 

which is, of course, the usual energy equation for laminar 
boundary-layer flow. In addition to equations (1) and (2), 
the equation expressing conservation of mass is needed: 

M + M =0 (3) 

dz by 

For a given gas the variation of /i and c v with temperature is 
known; hence, the foregoing system of three partial differen- 
tial equations is completed by the addition of the equation of 
state for a region of constant pressure 



ASSUMPTIONS 


2-percent error in the boundary-layer velocity profile; and 
hence, for all practical purposes, the mixing-layer velocity 
distribution calculations for Pr=l should be sufficiently 
accurate for air. 

Assumption 3 needs some explanation since the introduc- 
tion of a constant G differing from unity in the approximate 
relation between viscosity and temperature apparently has 
not been used in previous work. Usually G is taken as unity, 
and in such cases the approximation = {TjTJ) w gives 
reliable results for a fixed «, provided the free-stream tem- 
perature is restricted, to a certain range. By introducing the 
factor G, the approximating equation can be made to give 
the same value as a more exact equation at any deshed 



Figure I.— Schematic drawing of the flow. 


temperature in the mixing layer regardless of or w. 
Assuming that Sutherland’s equation 


ix f T Y* T a +S 
n„~VTj T+S 


(5) 


represents the true variation of viscosity with temperature, 
then the approximate equation can be made exact at any 
given temperature T by means of the relation 


C 



3/2 T„+S 

~T+& 


( 6 ) 


In order to solve the above system of equations, the follow- 
ing assumptions are made : 

1. c p = constant 

2. Pr=c v y/k— 1 

3. JP—c( -J-Yj where G is a constant depending on T 

v* x 1 »/ 

and T a 

The second of these assumptions if often made in calculating 
boundary-layer flows when only the velocity profile is desired 
and not the thermal characteristics. The difference between 
the boundary-layer velocity distributions for Pr = 1 and 
Pr= 0.73 is small, as is clearly shown by the numerical results 
of Emmons and Brainerd (reference 3). Since the mixing- 
layer and boundary-layer flows differ only in the boundary 
conditions and not in the differential equations, the effect of 
assumption 2 may be expected to be similar in the two types 
of viscous flow. At moderate supersonic Mach numbers, the 
use of Pr=l for air does not introduce more than 1- or 


In particular, if the approximating equation is linear in 
temperature («=1, thereby greatly simplifying the boundary- 
layer equations) and the viscosity is matched at the tem- 
perature T d , then the above equation gives 

°“Y T a % +8 

as the value of the constant G. By selecting C in this man- 
ner, rather than taking it as unity, a linear variation of 
viscosity with temperature then becomes an accurate ap- 
proximation in the inner part of the viscous layer, rather 
than in the outer part where the viscous stresses are less 
important. 

SOLUTION TO BASIC EQUATIONS 

As was first pointed out by Prandtl in reference 4, and 
later used to advantage- by Busemann and Crocco (refer- 
ences 5 and 6, respectively), the consequence of the assump- 
tion Pr= 1 when applied to boundary-layer flow is that the 
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temperature becomes a function only of the velocity. Hence 

c,T=m (8) 

By substituting this relation into equation (2) and using 
equation (1) in conjunction with the assumption Pr=l, it 
follows that the energy equation is automatically satisfied 
if the function/(u) satisfies the ordinary differential equation 

H? +i =° «•> 


Integrating this equation, using the boundary conditions, 


T= Ta> for u= U a 


gives 


T=T d for u =0 
f(u)=c P T=c p T i -~+~ 



( 10 ) 

( 11 ) 


It follows that 

fbu\ /d-iA /d«\ 

^ \bx)^ pl ' \by)~ pU 

and 

/du\ pu /bu\ 

•UA—p-UA 

Hence the momentum equation (1) in the (x, p) system 
becomes 


bu & / p bu\ 
p ° bx~bt V Up p» bp) 


( 13 ) 


This can be put in dimensionless form by introducing the 
variables 




as the relationship between velocity and temperature. 
Since the temperature determines the density, equation (11) 
also provides a means for calculating the density as a func- 
tion of the velocity. 

Following the method first given for incompressible flow 
by von Mises (reference 7 ) and later used for compressible 
flow by von Karman and Tsien (reference 8), a transforma- 
tion is made to a new set of independent variables (as, p), 
where p is the stream function. By using p as one of the 
independent variables, the continuity equation ( 3 ) is iden- 
tically satisfied, and the velocity components are given by 


u— 


Pm Z>P c== _P»d^ 
p by’ p bx 


( 12 ) 


Since the requirements of conservation of energy are ful- 
filled by equation (11), and conservation of mass by equa- 
tion (12), the only equation now remaining to be satisfied 
is the momentum equation (1). If a transformation were 
made to a completely new set of independent variables 
(s, p), the transformation formula would be 


P 


* 



P*= 


'P 

V vjjjoo 


( 14 ) 


Except for the parameter G appearing in the definition of p* 
and p*, these variables are the same as those used by Karman 
and Tsien (reference 8). Remembering that T 0 *, the 
free-stream total-temperature ratio, is given by 




MJ 


then the relation (equation (11)) between temperature and 
velocity can be written as 

T*= T a *-^- MJu**+ (T 0 *- T a *)u* ( 15 ) 

The momentum equation ( 13 ) becomes, using p*T*~ 1, 

*> u *- d del 

dr* bp* \ bp*) 


& bp b . bs d _pu d . bs b 
by~by bp~*by bs p a bp by ds 

d bp b_,bs ^ po_ , 5 s _d^ 

dr — dr bp^~bx bs~ p a bp'bx bs 


Setting s=r, 


bs . . ds , 

3— =0 and 5— =1 
by bx 


so that the transformation formulas are 2 


_d __ pu _d_ 
by~ Po bp 



3 The variables held constant in a differentiation process are explicitly Indicated in those 
cases where ambiguity could result if the subscript notation were not used. 


This is the basic equation which must be solved. The 
boundary conditions of the problem are such that no bound- 
ary layer exists at the point where mixing first begins. 
Under these conditions the velocity profiles will be similar 
at all points downstream of the origin, and hence the velocity 
u* will be a function only of some dimensionless variable f. 
This dimensionless variable must involve both p* and 2*, 
and must be zero at the origin of coordinates since the mixing- 
layer thickness is also zero at the origin. Therefore, let 

£=p* a x* i 

where a and 6 are pure numbers which must be determined 
by the condition that both sides of equation ( 16 ) for u* are 
functions only of the single variable f. Setting 

ff (f)su*T*“" 1 ( 17 ) 
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then the right side of equation (16) can be written as 


d / g du* df 5 


*) a**( 




)= 


from which it is obvious that in order for the right side of 
equation (16) to be a function only of {, it is necessary that 
a=l. With u*=u*(£) and a= 1, the partial differential 
equation (16) reduces to the ordinary differential equation 


bx 


du* 

W 


d / du*\ 
'~d£ \ 9 d£ ) 


Consequently, in order for the entire equation to be a func- 
tion only of f, it is also necessary that b=— b hence 


V = _ f 

■yfx* ^U a v a xG 


(18) 


and the ordinary differential equation for the velocity dis- 
tribution now reduces to 


£du* d / du*\ 

- -a?(* W 


2 d{ 


) 


(19) 


Equation (19) is the same differential equation that was 
obtained for boundary-layer flow in reference 8. It is a non- 
linear differential equation since g depends upon the velocity 
u*. This equation, however, can easily be converted into 
an integral equation which can be solved by the method of 
successive approximations. The conversion is made by 
temporarily assuming that fir is a known function of f (instead 
of u*) and formally applying the standard methods for 
solving first- order linear differential equations. The result is 


u 


where 


= <?i Pf df+O, 

Jo 9 

- pi# 

Jo 29 


F=e 


The boundary conditions are 

«*=1 at f = <» 
«*=0 at f= — oo 


( 20 ) 

( 21 ) 

(22a) 

(22b) 


Letting u 0 * be the value of u*. at f=0, equation (20) can be 
written as 

u*=Ci Ppdf+'Mo* (23) 

Jo y 

The constant Ci must satisfy two requirements 


1-V 


-C, 




f. i dt J~f* 


(24) 


Equation (23) is an integral equation for u*, since both F 
and g are functions of u*. By simply estimating a reasonable 


solution for u* as a function of f, a first approximation ]U* to 
the true solution can be calculated from 


lU 




* u F 


dt+ 0 Uo* 


The zero-order approximations 0 F and o 0 can be calculated 
directly from 0 u* by using equations (16), (17), and (21). 
If this process is repeated untill a given approximation is the 
same as the preceding one (to the degree of accuracy desired), 
and equation (24) is simultaneously satisfied, then the solu- 
tion to the problem is obtained. The iteration process turns 
out to b_e rapdily convergent, requiring two or three itera- 
tions to obtain the function t£*(f) accurate to within 1 per- 
cent, and about four or five iterations to obtain ix*(f) accu- 
rate to within a few tenths of 1 percent. 

In order to change the function «*(f) back to the physical 
coordinates ( x , y) a simple quadrature is necessary. By 
definition of the stream function, 

(I), *+($). i9+ <3£l * 


or 


1 dx+ dj . = pn dy _PZ dx 

Z \ X p co pee 

Hence, with a: held constant, integration gives 


(26) 


I U„ _ Ci T* 

V \v„xC Jo u* 

from which ufU* as a function of V^J~. can be determined 


It is to be noted that no graphical or numerical differentia- 
tions are needed at any point in the above iteration process, 
only quadratures are required. 

As is evident from equation (16), the assumption y*=CT*, 
that is, co=l, makes the momentum equation (16) independ- 
ent of temperature, and hence density. Consequently, with 
w=l, the solution to equation (16) in (x,\p) coordinates is 
independent of Mach number. For zero Mach number, 
T*= 1, (7=1, and 

(WSX..-I? 

Using this relation the solution in physical ( x,y ) coordinates 
for a linear variation of viscosity with temperature is obtained 
from the solution for zero Mach number by substituting 
equation (15) into equation (25). This gives 

^ u*d{+ (T 0 *- T/jf (26) 

where the integral in the second term on the right side is 
carried out forW*(f) corresponding to M.= 0. 
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RESULTS 

Numerical calculations of the Telocity distribution have 
been, made for the following cases : 

1 . p*=T*° r *; M a = 0 , 1 , 2, 3, and 5 

2. fi*=CT*; M a = 0, 1, 2, 3, and 5 

The various solutions for case 2 are obtained directly 
(equation (26)) from the solution for M a =Q without carry- 
ing out the laborious iteration process that is necessary to 
obtain solutions for case 1. All numerical results in physical 
coordinates (y-^UafvcX as independent variable) have been 
calculated for the case T d =T a . If the dead-air tempera- 
ture is radically different from the free-stream stagnation 
temperature, the proper velocity-distribution curves can be 
obtained b 3 r carrying out the integration indicated in equa- 
tion (25), since the function u*(f) in (x,\p) coordinates is 
independent of the thermal boundary conditions of the 
problem. 

Curves of u*(f) are shown in figure 2 for various Mach 
numbers. The corresponding curves in the physical plane 
are shown in figure 3 for the case fi*=T*°- 7s , and in figure 4 
for the case p*=CT*. In the latter two figures the familiar 
Blasius curve for the incompressible laminar boundary-layer 
flow is shown for purposes of comparison. The constant C 
that is used in figure 4 is determined by matching the viscos- 
ity coefficient at the temperature T S *=T 0 * according to 
equation (7). 

The particular curves shown in figure 4 apply for T a =400° 
R. Curves for any other temperature level T„ differ only 
in the constant factor C. 


CONCLUDING REMARKS 

A comparison is shown in figure 5 which illustrates the 
good agreement between velocity distributions calculated 
for the two approximations, m*= F* 0 - 76 and p*=CT*. At 
a Mach number of 2 or less, the curves for a linear variation 
of viscosity with temperature virtually coincide with the 
curves for a 0.76-power variation. For general use the 
linear approximation is reco mm ended since it gives results 
which are practically as accurate as the former, yet does not 
require a laborious iteration solution to be worked out for 
each Mach number. 

In general the laminar-mixing layer is several times 
thicker than the laminar boundary layer, as is illustrated in 
figures 3 and 4 where, for purposes of camparison, the Blasius 
profile is also shown. The rate of growth of mixing-layer 
thickness with increasing Mach number is somewhat larger 
than the corresponding rate of growth for a laminar boundary 
layer. The curves in figure 3 indicate a value^of roughly 

^-=1+0.UMJ 

o Jf=0 

for rate of growth of the mixing layer; whereas the corres- 
ponding value for a laminar boundary layer (see reference 
9, for example) is approximately 

J^- = 1+0.09.MJ 

°M=0 

This difference is to be expected since a larger percentage of 
low-density air exists in a mixing layer than in a boundary 
layer. 



Figobe 2.— Velocity distribution, in the (r, f) coordinate system. 
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Figure 3.— Velocity distribution for 0.76-power variation of viscosity with temperature. 


Figure 4. — Velocity distribution for linear variation of viscosity with temporaturo. 
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Figure 5.— Comparison, of velocity profile for a linear and a 0.76-power variation of viscosity 

with temperature. 


The foregoing statements, which indicate an increase in 
mixing-layer thickness with increasing Mach number, are 
based on the assumption that the Reynolds number (UaxjvJ) 
is held constant while the Mach number is varied. In most 
experimental apparatus the Reynolds number changes con- 
siderably with a variation in Mach number. Consequently, 
depending upon the particular experimental method em- 
ployed, the observed rate of mixing in the x direction may 
be either increased or decreased if the Mach number is 
increased. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif., Jan. 5, 1949. 
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